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ABSTRACT 

Context. Abundance anomalies observed in a fraction of A and B stars of both Pop I and II are apparently related to internal particle 
transport. 

Aims. Using available constraints from Sirius A, we wish to determine how well evolutionary models including atomic diffusion can 
explain observed abundance anomalies when either turbulence or mass loss is used as the main competitor to atomic diffusion. 
Methods. Complete stellar evolution models, including the effects of atomic diffusion and radiative accelerations, have been com- 
puted from the zero age main-sequence of 2.1 Mq stars for metallicities of Zq = 0.01 ±0.001 and shown to agree with the observed 
parameters of Sirius A. Surface abundances were predicted for three values of the mass loss rate and for four values of the mixed 
surface zone. 

Results. A mixed mass of ~ 10"' Mq or a mass loss rate of 10"" Mo/yr were determined through comparison with observations. Of 
the 17 abundances determined observationally which are included in our calculations, up to 15 can be predicted within 2 <t and 3 of 
the 4 determined upper limits are compatible. 

Conclusions. While the abundance anomalies can be reproduced slightly better using turbulence as the process competing with 
atomic diffusion, mass loss probably ought to be preferred since the mass loss rate required to fit abundance anomalies is compatible 
with the observationally determined rate. A mass loss rate within a factor of 2 of 10"''' Mg/yr is preferred. This restricts the range of 
the directly observed mass loss rate. 
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The brightest star in the sky, Sirius A, has all its main parameters 
such as mass, luminosity, age... relatively well determined, and 
its abundances hav e now been s tudied in great detail using Space 
Telescope data by iLandstreei (1201 lb . where more background 
information may also be found. 

In previous work, the then available observations of Sirius 
were compared with results from a grid of models calculated 
with atomic diffusi on and turbulence as a competing process by 
iRicher et al.l(l2000h . Given the range of abundances olDserved for 
any gi ven species, the ag reement seemed satisfactory. In a recent 
paper (IVick et al.ll20Tol) similar results were obtained through a 
similar approach but with mass loss as the competing process. 
Because the observational error bars were too large, these au- 
thors were unable to delineate which of mass loss or turbulence 
was responsible for its abundance anomalies. 

It is currently uncertain as to which process is most efficient 
in competing with atomic diffusion in A stars. For O and early B 
stars, mass loss is most likely the dominant process. It is clearly 
observed in those stars at a rate sufficient to obliterate the ef- 
fects of atomic diffusion. However in main-sequence A stars the 
expected mass loss rate due to radiative accelerations is smaller 
than in O stars by several orders of magnitude. Its presence i s 
likely only if it is started by another mechanism (lAbbottlll982h . 
and it might involve only metals (lBabel|[r995b . On the other 
hand, since surface convection zones are very thin, one expects 
little corona driven flux as observed on the Sun. It is then a pri- 
ori quite uncertain if A stars have any mass loss and the claimed 
mass loss rate for Sirius is an important observation. It is thus 



important to verify as precisely as possible if it is compatible 
with current observed abundance anomalies on the surface of 
Sirius. 

The measurement of the mass loss rate of Sirius is a diffi- 
cult observation and awaite d Hubble tel escope measurements of 
the Mgll resonance lines (Bertin et al. 1995). Their analysis of 
this spectral feature with a wind model leads to an uncertainty of 
0.5 dex on the mass loss rate if Mg is all once ionized. However 
there is an additional uncertainty related to the evaluation of the 
fraction of Mg that is once ionized. Their more credible evalua- 
tion of Mg ionization is based on setting the ionization rate equal 
to the recombination rate and leads tqj 



6 X lO"'"* < 



dM 
dt 



< 5 X lO-'^Mo/yr. 



(1) 



On the other hand, turbulence has often been used in stellar 
evolution calculations to explain observed abundance anomalies. 
In F stars of clusters, t urbulence could be responsible fo r the de- 
struction of surface Li jTalon & Charbonnelll998il2005h thereby 
leading to the so called Li gap. Turbulence could also be respon- 
sible for reducing abundanc e anomalies cause d by atomic dif- 
fusion on Am and Fm stars dTalon et alj|2006h . It can naturally 
explain the disappearance of abundance anomalies as rotation 
increases in those objects. It could also play a role for the Li 
abundance evolution in solar type stars (Pinsonneault et al. 198^ 



' This is slightly different from their Eq. [17] for which they had ne- 
glected the error bar given with their Eq. [7] but included an evaluation 
of the Mg II ionizati on based on a corrected LTE value coming from 
atmosphere models ISniiders & Lameril 19751) which do not seem ap- 
propriate for the wind region of interest here. 
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iProffitt & MichaudI Il99lh . It has however always been found 
necessary to use a number of adjustable parameters for its de- 
scription by physical models of turbulent transport and its role is 
uncertain. In this series of papers on the role of atomic diffusion 
in stellar evolution, turbulence was only introduced when mod- 
els with atomic diffusion led to anomalies larger than observed. 
Only one parameter was adjusted in order to control the influ- 
ence of turb ulence in limiting the size of anomalies: the surface 
mixed mass dRicher et al.ll2000tlMichaud et al.ll201 lb . 

It is possible to improve what we learn from the acccurate 
observations of Sirius by making more precise evolutionary cal- 
culations. Instead o f the g rid of solar metallicity models used in 
[Richer et all (l2000h and in lVick et all (12010), this paper uses two 
new series of models which were respectively computed with 
turbulence and with mass loss as the process competing with 
atomic diffusion. These two series are precisely converged to the 
known properties of Sirius, L, T^^fi, R, M and age. The origi- 
nal metallicity is determined and models with this metallcity are 
used to compare with observed abundances. This allows for a 
more precise and rigorous test than the calculations using grids 
of models. 

In this pap er stellar evol ution models with atomic dif fusion 
as described in [Richer et al.l (,2000) and in lVick et alj (1201 0) are 
used to determine the original metallicity of Sirius A using the 
age, radius and mass as constraints (Sect.|2]l. Using this metallic- 
ity and the determined parameters, complete evolutionary mod- 
els are then calculated (Sect.O, and the surface abundances are 
compared with Landstreet's recent observations (Sect.|4|. The 
results are discussed in Sect.jS] 

2. Original metallicity 

The fundamental parameters required to carry out stellar evolu- 
tion calculations, except for the original chemical composition, 
have been relatively well determined for Sirius A. The Hipparcos 
parallax can be used to determine its distance and when cou- 
pled with interferometry (Kervellaet al. 2003), to determine its 
radius, 1.711 +0.013 Rq. These authors also use the Hipparcos 
parallax to determine the luminosity from the magnitude and to 
refine the mass determination to 2.12 ±0.06 Mq. From the lu- 
minosity and radi us one c an obtain T^s - 9900+ 140 K while 
from spectroscopy jLe^i (1989) obtained T^ff = 9900+200 K 
The a ge of Sirius was discussed in Sect. 4. 1 of [Richer et al.l 
(120001) : using evolutionary time scales of Sirius B, they sug- 
gested 250±50Myr We adopt the slightly more re strictive 
range 225 ±25 Myr suggested by iKervella et al.l (l2003l) in their 
Sect. 2.2, where they argue that Sirius B would ha ve been more 
massiv e on the main-sequence than assumed by [Richer et al.l 
(l2000l) . Since the mass is well determined, we use a fixed mass 
of 2. 1 Mq, and then evolve models with a range o f metalhcities, 
using a scaled solar mix as given in Table 1 of iTurcotte et al.l 
(|1998|). Helium was adjus ted to the fitted solar value for model 
H of Turcotte et al ] (fT998l) as starting homogeneous composition 
for some of the calculations. For most of the calculations how- 
ever, the He mass fraction was corrected by AY = -0.02 dex 
because of the lower final metallicitjQ. 

In this paper a solar mass fraction, Xq, is often used in par- 
ticular to normalize both observ ed and calcu lated mass fractions. 
They are from Table 1 of iTurcotte et al.l (119981) and correspond to 



^ This 0.02 reductio n of Y for a 0.01 reductio n of Z is the same pro- 
portionality as used by IVandenBerg et al.ll2000l in building his Table 2. 
LiBeB were taken from the meteoretic values of iGrevesse & Sauvall 
119981 



the solar mass fractions at the birth of the Sun, more precisely the 
Zo(= 0.01999) for model H in Table 6 of that paper The surface 
solar abundances of metals today are some 10% smaller Those 
normalizing factors are different from the solar abundances used 
by Landstreet for comparative purposes. Since, in this paper, the 
same normalizing factors are used for both observed and calcu- 
lated quantities, they do not influence the comparison. 

Age, luminosity, Teff and radius are assumed well determined 
and are used as constraints to determine the original metallic- 
ity using models with turbulence. In Fig.[T] only models with 
Zo = 0.009, 0.010 and 0.011 are seen to satisfy all constraints 
within the predetermined error boxes. The radius is generally 
satisfied only for models younger than 250 Myr. Furthermore 
only models with the He mass fraction reduced by 0.02 satisfy 
the constraint on L. Models with higher Z and a solar He mass 
fraction do not satisfy the constraint from the radius as may be 
seen from the trend of the models with the solar He mass frac- 
tion. 

Below the deepest surface convection zone, the turbulent dif- 
fusion coefficient has been assumed to obey a simple algebraic 
dependence on density given, in the calculations with turbulence 
presented in this paper, by 



Dt = 10'*D(He)o'^° 



(2) 



where D(He)o is the atomic diffusion coefficient of He at some 
reference depth. Let AM = M, - M,- be the mass outside the 
sphere of radius r. For this paper, a series of models with differ- 
ent surface mixed masses (proportional to our parameter AMo) 
were calculated. More precisely calculations were carried out 
with 



PQ = p(AMo). 



(3) 



where Eq. (|3]l is given by the current stellar model. In words, 
in the calculations with turbulence reported in this paper, po of 
Eq. (|2]i is the density found at depth AM = AMq in the evolution- 
ary model. In practice the outer ~ 3 x AMq of the star is mixed 
by turbulence; for AMq = 10"^ '' M© the concentration of most 
species is constant for AM < lO"^'^ Mq. As one increases AMq, 
one defines a one parameter family of models. 

3. The models of Sirius A 

Two series of models were evolved for Sirius A. One in which 
th e process coinpeting with diffusion i s turb i ilence as described 
in [Richer et aP (120001) and Michaud etal.l (1201 ll) an d one in 
which it is mass loss as described in Vick et al.' ('2010'). In both, 
using opacity spectra from Iglesias & Rogers ( 1996), all aspects 
of atomic diffusion transport are treated in detail from first prin- 
ciples. 

The radiative accelerations of Mg, Si, Ca, Fe and Ni are 
shown in the upper row of Fig.|2| at ~232 Myr, the approximate 
age of Sirius A, in both a model calculated with mass loss (red 
curves) and one calculated with turbulence (blues curves). The 
two are practically superposed for log AM/M, > -5 but are sig- 
nificantly different for many species closer to the surface (that is 



3 The 
mula 

from the 

D(He) = 3.3 X lO-'^'r^ V[4pln(l + 1.125 x IQ-'TVp)] 



values of Z)(He)o actually used in this for- 
were always obtained — for programming con- 
venience — from the simple analytical approximation 

(in cgs 



units) for He in trace amount in an ionized hydrogen plasma. These can 
differ significantly from the more accurate values used elsewhere in the 
code. 



G. Michaud et al.: Sirius A: turbulence or mass loss? 



3 



10400 
1.50 



1.45 - 







1.40 



1.35 - 



10200 




9600 



9400 



^,^^r. 



solid curves From 200 to 250 My 

I ■ ■ ■ I ■ ■ ■ I ■ ■ ■ I ■ 



4.34 



4.32 - 



I 4.30 



00 

o 



4.28 - 



T — I — I — I — I — n — I — r 



0.011 -0.02 
0.011 -0.02 
0.010 -0.00 
0.010 -0.02 
0.010 -0.02 





9600 - 



Age (Myr) 

Fig. 1. An HR diagram and the time evolution of Tetf, L, logg, R and Zjurf are shown. Observationally determined +lcr intervals are 
shown for L, Teff, g and R. For Teft the spectroscopically determined error bar is in black while that determined using luminosity and 
radius is in red (see the text). Each model is color coded and identified on the figure. The adopted acceptable age range is from 200 
to 250 Myr (see the text). In the HR diagram, the part of the curves between 200 and 250 Myr is solid; it is dotted outside of that 
interval. All models were calculated with turbulence. Models with mass loss could not be distinguished, in the HR diagram, from 
those of the same original mass and composition calculated with turbulence. A subset of the curves is shown on Fig. lA.ll 



log AM/M» < -5). In the lower row the corresponding internal 
concentrations for the model with mass loss and for that with 
turbulence are shown. When there is a difference between the 
brad's for the two models they are caused by effects of saturation, 
as may be seen by comparing the abundances in the lower row. 
The X(Fe) and X(Ni) are larger in the model with turbulence for 
log AM/Mj < -5 than in the model with mass loss; the reduction 
of the photon flux at the wavelengths where Fe and Ni absorb the 
most is by a larger factor when the abundance is larger and so the 
brad's are smaller in the model with turbulence. The large under- 

abundance of Ca at log AM/M» 7.5 in the wind model causes 

the larger grad(Ca) there, but the large underabundance is also 
caused by the maximum of ^rad(Ca) as discussed in Sec. 5.1.1 of 
fVick et al., (.2010.) to which the reader is referred for a detailed 



discussion of the interior wind solution. While the surface abun- 
dances of say Fe and Ni in the model with mass loss are within 
0.3 dex of those in the model with turbulence, their interior mass 
fractions differ by a factor of about 5 for -7 < log AM/M, < -5. 

Figures lBTSl and lBTTI contain results for all calculated species; 
they are are shown in the online AppendixlBl 

4. Surface abundances 

Given the constraints imposed by age, Teff, L and R, and 
that the mass is 2.1 Mq, the original metallicity is fixed to 
Zo - 0.010 + 0.001 (see Sect.©. Specificafly, the luminosity 
(middle left hand panel of Fig.© determines Zq = 0.010 + 0.001 
and then the radius determines the acceptable age to lie between 
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Fig. 2. Upper row Radiative accelerations for five atomic species, after 232 Myr of evolution in a model with a mass loss rate of 
10"'^ Mg/yT (red curves) and in a model with turbulence (blue curves). The dotted lines represent gravity. Lower row Corresponding 
mass fraction in the model with mass loss and in that with turbulence. The dashed lines are the original values. They correspond to 
an original metallicity of Zq = 0.01. 
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Fig. 3. Observed and predicted abundances on Sirius A as a function of Z, the atomic number, left in the model with turbulence 
calculated with an original metallicity of 0.009 for three slightly different values of turbulence and right in the model with mass loss 
with an original metallicity of 0.010 calculated with four different mass loss rates. The model name Z0.009_dY-0.02_mtbl.O stands 
for a model calculated with Zo = 0.009, AY ^ -0.02, and AMq = 1.0 x lO"*" Mg. The model name MassLossJZ0.010W2E- 13 stands 
for a model calculated with dM/dt - -2x lO^'^* Mg/yi. All dotted lines represent models of about 200 Myr. 
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Fig. 4. Observed and predicted abundances on Sirius A as a function of Z, the atomic number, in the model with turbulence for 
two slightly different values of original metallicity, 0.010 and 0.011, and in the model with a mass loss rate of 10''-^ Mo/yr and 
a me tallicity of 0.01 (see Fig.[3]for model labeling definitions). In the left panel, the two sets of observations are from lLandstreetl 
12011 : the black circles are from his observations while the pink triangles are an average over all recent observations (se e the text). In 
the right panel , the averaged valu es are replaced by the act ual data points o f each observer as given in Table I ofiLandstreet (1201 Ih 
where c ircles. iLandstree t[ (l2011h : in verted open fnawg / e^. iLambert et akl (Il982|): inverted three-po int st ars. iLemkd (ll99Cj): blu e 
squares, lOiu et all (l20'qi|i diamonds, jHill&Landstreef ("iggSb; asterisks, Hui Bon Hoa et al h99lh; pfa5. lRentzsch-Holml (Il99"7h : 
upright open fn'gngfe^. lHolweger & S turenburg ( 1993); pink squares. .Si?idaknne & Ueta (1989) . See the text for the explanation of 
the differences between the error bars in the right and left panels. 



200 and < 230 Myr (middle right hand panel). There only re- 
mains mass loss rates or mass mixed by turbulence that may be 
varied to define a range of predicted abundances that can then be 
compared with observations. 

Evolutionary models were calculated for AMq - 1 .0, 1 .4 and 
2.1x10"^Mq (see Eq. Q) and for mass loss rates of 0.5, 0.7, 
1.0 and 2.0 xlO"'"* Mo/yr. On Fig. [3] predictions from some of 
them are compared with observations from column 2 of Table 1 
of ILandstreet (2011). On Fig |4] data from the other columns of 
his Table 1 are also used in order to present a picture of the un- 
certainties of the observations, as briefly discussed below. 

In the left panel of Fig. [3] results are shown for the case 
Zo = 0.009 at 200 Myr (dotted line segments) and at 220 Myr 
(solid line segments) which is the age interval over which all 
constraints are satisfied according to Fig.[T] In practice, even 
though shown for all cases, the dotted and solid segments are 
barely distinguishable and merely widen the dots. Note the as- 
sumed original composition at age 0.0 Myr in light blue on each 
panel of Figs.[3]and|4] As the mass mixed by turbulence is de- 
creased from 2.1 to LOxlO'^M©, the surface abundance of el- 
ements supported by grad (e-g- most Fe peak elements) and the 
underabundance factor of sinking elements both increase. In the 
left panel of Fig.|4]similar results are shown for original metallic- 
ities of Zo = 0.010 and 0.01 1 with AMq = 1.4 x 10"^ Mq. Within 
the range of original metallicities acceptable according to Sect.|2] 
(from Zo = 0.009 to 0.01 1), the level of agreement between pre- 
dicted surface abundances and observed ones does not change 
much although the Zq = 0.010 case is slightly favored. Given 
the number of observed species, there is in practice little room 



for adjustment: as one may see from the left panel of Fig. [3] the 
Fe peak abundances favor the lower value of the mixed mass but 
the abundances of He, O, S and Ca rather favor the larger value. 

Predictions for four mass loss rates and a metallicity of 
Zq - 0.010 are shown in the right panel of Fig. [3] Only one 
metallicity is shown since the effects of changing metallicity 
over the acceptable range of Zo = 0.009 to 0.011 are small 
as briefly mentioned above for the turbulence models. The 
iron peak elements would probably favor a mass loss rate of 
0.5 xlO~'^ M^/yi, but the lower mass species then disagree 
more strongly, and the better compromise appears to be the 1 .0 
xlO"'^ Mo/yr case. The iron peak elements show approximately 
the same sensitivity to the mass loss range as to the mixed mass 
range illustrated in Fig.[3]but He, O and S are more sensitive to 
the mass loss ratfl 

In the left and right panels of Fig.|4]are shown the same three 
theoretical models (two calculated with turbulence and one with 
mass loss) com pared to two differe nt presentations of the results 
from Table 1 of ^L andstreetl (1201 li") . On the left panel are shown 
both his determinations (his column labeled Li 7) and, as sepa- 
rate data points, his determinations averaged with those of the 
other observers he lists in his table, except for a few values he 



* As the mass loss rate is reduced, the settling velocity becomes 
closer to the wind velocity in magnitude. This tends to amplify differ- 
ences in settling velocities caused, among other factors, by small mass 
differences. For instance the abundance variation of "^He is a factor of 
1.8 larger than that of ^He for the 0.5 xlO"'^ M^/yr case but only a fac- 
tor of 1.3 larger in the 10"'^^ Afo/yr case (see the right panel of Fig.|3j. 
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Fig. 5. Color-coded interior concentrations for the same Zo - 0.010 models as in Fig.|4]at ~ 233 Myr left in the model with 
turbulence and right in the model with mass loss. The radial coordinate is the radius and its scale is linear, but the logarithmic value 
of the mass coordinate above a number of points, log AM/M», is shown on the left of the horizontal black line. The concentration 
scale is given in the right insert. Small circles near the center of both models mark the central convection zone. While the surface 
abundances are very similar, as seen in Fig.|4] the interior concentrations are quite different between log AM/M, - -3 and the 
surface. 



argues are erroneous or too uncertain to be worth including (his 
column labeled mean). The error bars are standard deviations of 
the observations listed in his table. These do not include con- 
tributions from the error bars of the various authors. The actual 
values of the different observers listed in his Table 1 are shown in 
the right panel of Fig.|4] The values which Landstreet excluded 
from his averages are not shown but all others are shown. One 
can argue that the difference in abundance values between the 
various observers is a better estimate of the uncertainty of the 
abundance determinations than the mean value with associated 
standard deviation shown in the left panel. Following a discus- 
sion with Landstreet (private communication) the error bars for 
his points (from his column Lll) were slightly increased for 
the right panel of Fig.|4]only. So all cases where, in the "Lll" 
column of Table 1, he had 0.1, we increased to 0.15 because 
this is closer to the actual dispersion found for dominant ions 
with many Unes. Noting that there is some additional uncer- 
tainty (about 0. 1 dex) due to imprecise fundamental parameters 
and microturbulence, another 0.1 dex was added in quadrature to 
all sigmas, giving a minimum sigma of 0. 18. 

The right panel of Fig.|4]gives our comparison between the- 
ory and observations. The agreement is not perfect for any value 
of turbulence or mass loss. In fact the difference between the 
mass loss and turbulent models is very small for atomic species 
lighter than Cr and one may question if abundances alone can 
really distinguish between the two given the error bars. Of the 
species included in our calculations, abundances were deter- 
mined observationally for 17 atomic species and upper limits 
for 4. For the model Zq = 0.010 with AMq = 1.4 x IQ-^ Mq 
one counts 8 species within 1 cr and an additional 6 within 2 
cr. These numbers are respectively 8 and 7 for the Zq - 0.011 
model with the same turbulence. One counts respectively 8 and 
5 for the Zo = 0.010 mass loss model with a mass loss rate of 
1.0 xlO"'"' Mo/yr. In addition, of the 4 upper limits, 3 are com- 
patible with predictions. While the trend is right, the agreement 
is not so good for Ti, Cr and Mn. 



The interior concentrations in the two Zq = 0.01 models 
of Fig. |4] are shown in Fig.|5] While the surface abundances are 
quite similar in the two, the interior concentrations are quite dif- 
ferent for log AM/Mt < -3 (see Sect.|3]l. 

If one compares these results with Fig. 20 of IVick et al] 

(1201 01) . one notes that the same mass loss rate of 
1.0 xlO~'^ Mo/yr had been found to lead to the prediction 
closest to observed abundances. While the age assumed for 
Sirius is about the same in the two p apers, the larger mass of 
the models used in IVick et al.1 (12010) causes them to be more 
evolved and, so, have a smaller gravity at a similar age. This is 
an important difference, when one compares two models with 
approximately the same Teff and thus the same giad's. Another 
difference comes from the original Zq, which is 0.010 for the 
mass loss models of Fig.|4]of this paper, but 0.02 for those of 
Fig. 20 of lVick et al.1 (I201Q) . 



5. Conclusion 

Using observationally determined stellar parameters for Sirius A 
one first fixed a metallicity and He mass fraction (Sect|2]i. Then, 
expected surface abundances were predicted as a function of ei- 
ther a surface mass mixed by turbulence or of a mass loss rate. 
Of the 17 abundances determined observationally, up to 15 can 
be predicted within 2 cr, and 3 of the 4 determined upper limits 
are satisfied. The three atomic species, B, N, and Na show the 
strongest disagreement. 

While the origin of the assumed turbulence is not de- 
termined, it could be either shear induced by differential 
rotation ( Talon et al. 2006) or it could be gravity waves 
(iTalon & Charbonnellll998l) . If the mam competmg process is 
mass loss howe ver, it has the gr eat advantage of being al- 
ready observed dBertin et al.lll995l) . The mass loss rate of 1.0 
xlO"'-' Mo/yr that best reproduces abundance observations is 
slightly larger than the lower limit of 6 x 1 M^lyx determined 
from asymetries of Mg II lines using ST observations by these 
authors (see Eq. IT]). Their estimate based on corrected LTE, 
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their Eq. [12], however gives a mass loss rate between 5.0 xlO" 
and 5.0 x 10"" Mo/yr which would lead to practically no surface 
abundance variations during evolution in contradiction with the 
obse rved Sirius abundances (see Figs. 11 and 20 of I Vick et al.l 
12010. for a calculation with 1.0 xlO"'^ Mg/yr). Our results then 
support the arguments presented above in Sect.[T]in favor of their 
estimate based on radiative ionization fraction. While by them- 
selves, abundances do not favor mass loss or turbulence as the 
competing process, the agreement with the observationaUy de- 
termined mass loss rate favors mass loss. 

While there probably still remains some uncertainties in the 
observationaUy determined abundances, the disagreement be- 
tween observations and calculations points to some weaknesses 
of the models. The first one may be related to the presence of 
Sirius B. Even though it is quite a wide pair, it had been sug- 
gested by Richer etal. (2000) that the disagreement with the 
C and N observations could be caused by the transfer of ma- 
terial from the fo rmer primary. This was further discussed by 
iLandstreefl (1201 ih who also suggested that it could simultane- 
ously explain the difficulty with the B upper limit. Sodium could 
also be affected, jus t as it is affected in the globular cluster M4 
(iJVlarino et alJl201 ih . However the extent of the mass transfer in 
Sirius remains uncertain. 

Even if it is tempting to accept mass loss as the most 
important mechanism competing with diffusion in slowly ro- 
tating stars, it is also observed that abundance anomalies are 
much less important in rapidly rotating stars. An other mech- 
anism linked to rotation must then be involved. Either rota- 
tion dr iven turbulence ( Talon et al. 2006) or meridional circu- 
lation dCharbonneau & Michaudiil99li) could progressively re- 
duce abundance anomalies as rotation increases. The weak de- 
pendence of abundance anomalies on the rotation rate could be 
due to its effect becoming larger than those of mass loss only as 
one approaches the 100 km/s limit of v sin / (Abt 2000) for the 
Am star phenomenom. 

In relation to Fig.|5j it was suggested in Sect.|4]that asterosis- 
mology tests could distinguish between mass loss and turbulence 
as the competing mechanism for slowly rotating stars. While de- 
tailed evolutionary calculations with meridional circulation have 
not yet been caiTied out, one expects that laminar meridional cir- 
culation would lead to internal metal distributions, in 3 dimen- 
sions, similar to those that mass loss leads to, in 1 dimension, 
since both are advective and not diffusive processes. This opens 
the possibility of distinguishing between meridional circulation 
and rotation induced turbulence using asterosismology. 

For simplicity, these calculations were carried out with un- 
differentiated mass loss throughout Sirius A's evolution as seems 
appropriate for Am stars. However a 2. 1 Mq star starts its main- 
sequence at Teff ~ 10500K (see Fig.[TJ where no H convection 
zone is present and which is probably within the HgMn domain. 
What would have been the mass loss rate then? The observa- 
tion of Hg isotope anomalies on Hg Mn stars suggests that the 
mass loss would be differentiated ( see iMichaud et al.ll 19741 and 
Sect. 4 of IMichaud & Richedl200 8^. How would this affect sur- 
face abundances during later evolutionary stages such as reached 
by Sirius A? 
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Appendix A: Properties of models 

Appendix B: Radiative accelerations and interior 
concentrations for all species 
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Fig. A.l. An HR diagram and the time evolution of Tetf, L, \ogg, R and Zsm-t are shown. Observationally determined +lcr intervals 
are shown for L, Teff, g and R. For Teff the spectroscopically determined error bar is in black while that determined using luminosity 
and radius is in red (see the text). Each model is color coded and identified on the figure. The adopted acceptable age range is from 
200 to 250 Myr (see the text). In the HR diagram, only the part of the curves between 200 and 250 Myr is shown; see Fig.[T]for a 
more complete figure. All models were calculated with turbulence. 




log (AM/M.) 

Fig. B.l. Radiative accelerations after 232 Myr evolution in a model with mass loss of 10 Mo/yr (red curves) and in a model with 
turbulence with AMq = 1.4 x 10"^ (blue curves) for all calculated atomic species. The sharp mimima in the ^i-ad(Ti), grad(Cr), 

and grad(Mn) curves at log AM/M» 10 appear surprising but have been verified to be a systematic property of the OPAL atomic 

data. 
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Fig.B.2. Mass fractions in the models with turbulence (blue curves) and mass loss (red curves). See the caption of Fig. lBTI 



